Nanoscale electronic order in iron pnictides 
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The charge distribution in ^FeAsOi-ajFj, (/?=La, Sm) iron pnictides is probed using As nuclear quadrupole 
resonance. Whereas undoped and optimally-doped or overdoped compounds feature a single charge environ- 
ment, two charge environments are detected in the underdoped region. Spin-lattice relaxation measurements 
show their coexistence at the nanoscale. Together with the quantitative variations of the spectra with doping, 
they point to a local electronic order in the iron layers, where low- and high-doping-like regions would coexist. 
Implications for the interplay of static magnetism and superconductivity are discussed. 

PACS numbers: 74.70.Xa, 75.25.Dk, 76.60.-k 



Strong electronic coiTelations, notably present in several 
transition metal oxides, give rise to a broad range of exotic 
electronic states, in particular high-temperature superconduc- 
tivity. Major interest has followed from the recent discovery 
ilt] of the latter in iron pnictides, with critical temperatures 
second only to those of copper-based materials. Iron pnic- 
tides feature electronically active FeAs layers with a multior- 
bital character, separated by RiOi^^F^) layers (7?=rare earth) 
for the "1 1 1 1" family studied here. On doping with electrons 
through fluorine substitution, the antifeiTomagnetic order of 
the parent material RFeAsO is suppressed and a superconduc- 
tivity region shows up in the phase diagram. While the itiner- 
ant nature of the undoped compound contrasts with the Mott 
insulating character of the undoped cuprates, the renewed 
proximity of static magnetism and superconductivity is fuel- 
ing the idea that competition with another ground state is cru- 
cial to high-temperature superconductivity. However, reports 
in underdoped pnictides are so far conflicting, with micro- 
scopic i^-B 1 or mesoscopic flBI coexistence, a second-order 
boundary |6|], or a first-order boundary |7]. While this can be 
accounted for to some extent by differing sensitivities of the 
experimental probes, in particular to disordered magnetism, it 
remains unclear whether intrinsic electronic inhomogeneities 
and an associated order, short-range or more, can show up 
such as in cuprates |8- 10], nickelates |fTll[l2^] and manganates 
|13lll4|]. In the presence of such inhomogeneities, both the 
magnetic ordering and superconductivity ground states may 
be altered, impacting their interplay or competition. 

Any intrinsic inhomogeneities may be sensitive to the pres- 
ence of impurities within the FeAs planes, as in Co-doped 
BaFe2As2, or to the electrostatic potential associated to the 
dopant layers. As the 1111 pnictides correspond to the case 
where the dopants are the farthest away from the FeAs layers 
and with rare earth layers in between, they represent the situa- 
tion of minimal influence of the dopant layers on the intrinsic 
electronic properties. Therefore, we investigated the charge 
distribution in /?FeAs(Oi_a;Fj;) {R=La, Sm) 1111 compounds, 
over the 0<a;<0.15 doping range, using nuclear quadrupole 
resonance (NQR) measurements. Our study shows the pres- 
ence of an electronic inhomogeneity in the underdoped region 
of the phase diagram, ascribed to a nanoscale electronic order 
and independent from the low-temperature behavior. Reveal- 



ing that the phase diagram has to be understood beyond a ho- 
mogeneous picture with ground state competition, this result 
also supports the idea that local electronic order is a common 
feature of strongly-correlated systems, notably in transition 
metal oxides jlSll . 

All NQR measurements were carried on powder samples, 
which were prepared and studied using x-ray, susceptibility, 
resistivity, and /iSR measurements as described previously 
llTlfieHlsll. Undoped La and undoped/4% (a;=0.04) Sm sam- 
ples were shown to display magnetic ordering at T'7\/=138 K 
(La) and Tm= 137/95 K (Sm), with the accompanying struc- 
tural transition at T5=156 K (La) and Ts=160/140 K (Sm). 
The 5%/7.5%/10%/15% (La) and 6%/8%/10% (Sm) samples 
display bulk superconductivity below rc=20/21/26/l 1 K (La) 
and rc=36/45/52 K (Sm). To probe the atomic populations 
coiTesponding to potential electronic inhomogeneities, NQR 
takes advantage from the fact that a nucleus with a nuclear 
spin />l/2 features a finite electric quadrupole moment. In 
the presence of a finite electric field gradient (EFG) at the 
nuclear site, the degeneracy of the corresponding nuclear en- 
ergy levels is lifted, and probing by radiofrequency iiTadiation 
can be performed. As the EFG stems from the surrounding 
charge distribution, peculiarities of the latter can be inferred 
from the determination of the EFG histogram in the sample. 
Since 7=1/2 for iron, the ^^As nuclei (I —3/2) were used as 
NQR probes. One measures then the quadrupole frequency 
vq oc Vzz + where Vzz and 77 are respectively the 

highest eigenvalue and the asymmetry of the EFG tensor The 
proximity of the As ions to the iron layers helps to retain high 
sensitivity to electronic changes, which may be further helped 
by their large polarizability 1 19]. 

We start by examining the static properties, i.e., the his- 
togram of the electric quadrupole frequency vq. A standard 
■|-r-7r spin echo sequence was used, with t=24-34 /is and 
integration of the full echo. The r values, as well as the pulse 
repetition rate, were checked to be small enough that spin-spin 
or spin-lattice relaxation contrast is limited to at most a few 
percent of the relative line intensities. Spectra were connected 
for frequency-dependence of the intensity. On doping (see 
Fig. [TJ, the frequency distribution shifts to higher values, in 
keeping with previous observations l20l l2l[l and in disagree- 
ment with local-density approximation calculations 
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Figure 1: (Color online) ''^As NQR spectra of RFeAs(Oi-:rF^) 
(R=La, Sm). "OPT" and "OVD" refer to optimally doped and over- 
doped samples. Fits including up to three (La) or four (Sm) Gaus- 
sians are shown as full lines, with the two-Gaussian fit for a;=0.05 
(La) detailed as an example. 
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Figure 2: (Color online) (Upper panel) Doping dependence of the 
spectral peak frequencies, taken at the centers of gravity of the low- 
and high-frequency peaks of Fig. [T] For each rare earth, only points 
taken at similar temperatures are shown. Dotted lines are visual 
guides. (Lower panel) Doping dependence of the relative weight 
of the high-frequency peak. The solid blue line is a linear fit for 
0.04<a;<0.08. Dotted lines are visual guides. 



In the undoped limit, the single narrow line agrees with a sin- 
gle well-defined charge environment for all As nuclei. The 
line is broadened in the optimally-doped or overdoped limit 
(La 10%/15%, Sm 10%), likely reflecting the structural disor- 
der of the fluorine dopants and some moderate fluorine con- 
centration inhomogeneities. In the underdoped region (La 
5%/7.5%, Sm 4%/6%/8%), two fairly broad peaks are ob- 
served, with further structuring of the high-frequency peak for 
Sm samples. Spectral weight is transferred directly to high 
frequency on doping, with a roughly linear dependence ex- 
trapolating to the optimal doping (see lower panel of Fig.|2]l. 
Crucial questions are then the difference in nature and the 
scale of coexistence of the two corresponding sets of charge 
environments. 

A direct explanation would be phase separation on a macro- 
scopic or mesoscopic scale, with the difference in peak posi- 
tions (see upper panel of Fig. |2]i indicating low and high dop- 
ing regions. Beyond incompatibility with initial x-ray charac- 
terization, this can be tested using Tj^^ spin-lattice relaxation 
rate measurements in the La samples. For low-energy spin 
fluctuations over all wave vectors q. 
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where "^^Ahf is the As hyperfine form factor, x" is the dy- 
namic spin susceptibility, and w, / is the irradiation frequency. 



The relaxation was measured using the inversion recovery se- 
quence, fitting the recovery curve to 

M{t) = Mo(l - /exp(-(3i/Ti)^)), 

where A/q is the magnetization at thermal equilibrium, / ac- 
counts for incomplete inversion, and A accounts for spreading 
of the Ti value. If set free, A tended to drift below unity be- 
low Tc, with no significant effect on the extracted Ti value. 
Figure [3] shows for each spectral peak the temperature depen- 
dence of (TiT)^^ (upper panel) and of the same T^^ data, 
rescaled about Tc (lower panel). While {TiT)^^ in the un- 
doped material tends to diverge on approaching the magnetic 
transition as expected, all other peaks show no signature of 
magnetic ordering. They reflect however a superconducting 
transition at low temperature as seen from the rapid decrease 
of the relaxation, with a ^ behavior broadly consistent with 
observed power laws l23ll24ll . For the underdoped samples in 
the paramagnetic state, the two spectral peaks feature similar 
(TiT)^^ behaviors, very different from the progressive sup- 
pression of low-energy excitations observed at optimal dop- 
ing. If phase separation would occur even on a rather small 
scale (several nanometers or more), the volume fraction cor- 
responding to each peak would exhibit fluctuations specific to 
the then necessarily different doping levels. In light of the 
spectra, the relaxation contrast should then be much larger, 
with the high-frequency peak relaxation closer to developing 



Pubfished version: Phys. Rev. Lett. 104, 097001 (2010) 



3 



a decrease similar to that at optimal doping. Here, the sim- 
ilar weak Curie- Weiss behavior above Tc and the moderate 
difference in amplitude show a sharing of electronic proper- 
ties, i.e., the coexistence of the two charge environments at the 
nanoscale. The adequacy of using the same Tc to rescale the 
Tj^^ data for the two peaks at a given doping (lower panel of 
Fig-E) is consistent with this conclusion, at least for a;=0.075. 
Note that a comparable study cannot be performed on the Sm 
samples as Tj^^ is swamped by the contribution of the Sm 
magnetic moments. However, beyond the similarity to the 
La spectra, the increase of both peak frequencies on trans- 
ferring spectral weight (see upper panel of Fig. ^ is clearly 
inconsistent with a mere change in the proportion of two well- 
separated phases, and suggests mutual influence at the local 
level. 
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Figure 3: (Color online) ""As spin-lattice relaxation in 
LaFeAs(Oi_a;Fi:). (Upper panel) Temperature dependence of 
(TiT)-^ with T-^ the spin-lattice relaxation rate as measured on 
each peak of Fig. [T] with closed (open) symbols corresponding 
to low (high) frequencies. The measurements are done at the 
T-dependent peak frequency (undoped), 9.7 and 10.6 MHz (both 
5% and 7.5%), and 11 MHz (10%). It was checked on the 5% 
sample that the spectrum is preserved on crossing Tc- (Lower 
panel) Temperature-dependence of T^^ at low temperature, with 
horizontal scaling by Tc (as determined from initial characterization) 
and vertical scaling to obtain coincidence about Tc- Symbols for 
the low-frequency peak in the underdoped samples are oversized for 
emphasis. Dashed lines are visual guides. "SC" and "PM" refer to 
superconductivity and paramagnetism. 



Crucial questions are the origin of this coexistence and the 



spatial charge arrangement corresponding to it. A possibil- 
ity is fluorine control of the two charge environments at the 
As sites through direct electrostatic effect or local doping of 
the FeAs layers. To properly account for the weights of the 
peaks, each fluorine must influence on average roughly nine 
As ions (see lower panel of Fig. |2]i. This rules out a one-to- 
one effect of each dopant on the nearest arsenic, contrary ll25ll 
to the effect of Sr^"^ on Cu ions in La2-a;Sr2:Cu04. If a single 
fluorine would affect a patch of As ions, the linear growth of 
the high-frequency weight with doping suggests that no siz- 
able overlap of these patches develops as the doping rises. In 
particular, the fact that the linear behavior extrapolates well to 
the optimal doping indicates that in the latter case there would 
be a perfect arrangement of the patches, i.e., fluorine order- 
ing. While dopant ordering is known to occur in compounds 
such as sodium cobaltates |26], there is no experimental hint 
of it in pnictides. Therefore, the inhomogeneity likely arises 
not from the dopants but from an electronic instability in the 
FeAs layers, due to competing interactions. 

The resolved spectral features and their smooth doping de- 
pendence indicate rather well-defined and reproducible lo- 
cal electronic environments, static on the slow (microsecond) 
time scale of NQR. As seen on the upper panel of Fig. |2] 
the low-frequency and high-frequency environments are rea- 
sonably well connected, respectively, to the undoped and 
optimally-doped or overdoped environments. Together with 
the linear spectral weight transfer extrapolating to the optimal 
doping, this suggests a local electronic order with a varying 
ratio of low-doping-like and high-doping-like regions. Ex- 
pected to feature larger intrinsic magnetic fluctuations, the 
low-doping-like regions would account for the relaxation re- 
sponse of the whole system. Since doping corresponds to 
changes in the iron 3d orbital occupancies, it must, however, 
be noted that what appears to be a difference in total occu- 
pancy (charge ordering) around As sites may also have an or- 
bital character (orbital ordering) 13311 . This scenario would 



be of interest in light of the argued link 11271 12811 between or- 
bital ordering and static magnetism. More generally, the pres- 
ence of a local electronic order is reminiscent of the situa- 
tion observed for instance in cuprates and manganates 
lITii [14I1 . such as stripe or checkerboard order, static or dy- 
namic. This supports the widespread presence of electronic 
inhomogeneities in correlated systems llTsll . even in the pres- 
ence of a homogeneous ground state. 

Being common to all underdoped samples, the local order is 
no direct explanation for any ground state coexistence. How- 
ever, our data show that Tc remains high despite significant 
magnetic fluctuations, just as the superfluid density is com- 
paratively more depressed [7], possibly indicating a proximity 
effect in nonintrinsically superconducting low-doping sample 
regions. On the contrary, on doping, the static magnetism 
tends to be disordered QII^Ij with a large reduction of Tm- 
Therefore we propose that, beyond direct ground state compe- 
tition, the electronic inhomogeneity of the underdoped region 
of the phase diagram influences the transition from static mag- 
netism to superconductivity, as shown schematically on Fig.|4] 
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Figure 4: (Color online) Electronic phase diagram for the 1111 pnic- 
tides. TmITc denote the magnetic and superconducting transition 
temperatures. A local electronic order develops in the underdoped 
region, with superconductivity being favored over static magnetism. 
At the top are shown typical NQR spectra, taken from Fig.[T] 



Starting from high doping, superconductivity is unhindered if 
not helped by the local order and disappears [2,, ^ JJ only at 
its low-doping end, where the high proportion of low-doping- 
like regions would allow static magnetism to shoot up, before 
recovering electronic homogeneity close to the undoped limit. 
Note that, as is the case in other systems llsoll . the local or- 
der occurs where a quantum critical point could have been 
expected in the phase diagram. A crucial question is then 
whether the local environments correspond to specific points 
of the phase diagram, as is seen fsT] in La2-a;Sra;Cu04+j, 
with separation in stripe-ordered and optimally doped super- 
conducting regions. 

In conclusion, using As NQR, we show the presence of two 
charge environments at the nanoscale in underdoped 1111 iron 
pnictides. These can be ascribed to a local electronic order 
in the iron layers, in line with electronic inhomogeneities ob- 
served in other correlated compounds. We propose a phase di- 
agram in which this local order serves as the electronic back- 
ground for the interplay of static magnetism and supercon- 
ductivity, beyond direct competition. An intriguing scenario, 
to be addressed in future work, is the possible orbital nature 
of this local order, which could suggest that a primary factor 
in the phase diagram is real-space orbital competition, where 
sensitivity to structural details would yield seemingly differ- 
ent phase diagrams. In light of the reported importance of the 
Fermi surface topology [321 . this would represent a significant 
change of perspective on this system. 
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